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Abstract — In this research, analytical method of water 
wave dynamics was developed using wave constant Gat 
velocity potential of Laplace equation solution result, 
where at small amplitude theory the wave constant was 
eliminated with wave amplitude. The assumption of 
sinusoidal wave is maintained. 

The calculation methods that were formulated include 
wave number calculation method and shoaling and 
breaking modeling. Using the wave constant, then 
equilibrium equation will be met more accurately. 

Keywords — Wave constant G, sinusoidal water wave 
surface, wavelength, shoaling-breaking. 

I. INTRODUCTION 

Small amplitude wave theory (Dean (1991)) was 
formulated with an assumption that wave amplitude is 
very small so that the elevation of water wave surface at 
the execution of Bernoulli equations at the surface is 
considered to coincide with the still water level elevation 
(zero elevation). Using this method, relation equation 
between wave constant G and wave amplitude A was 
formulated. Afterwards, wave constant G at velocity 
potential equation was eliminated with wave amplitude A. 
Bearing in mind that there is no wave amplitude at the 
velocity potential equation, then the accuracy of wave 
analysis by eliminating wave constant G requires accurate 
relation equation. Particle velocity equation which is a 
differential of velocity potential equation is a function of 
wave constant G where the particle velocity is used in 
various calculations and at equilibrium equation. Hence, 
eliminating wave constant G with wave amplitude 
requires an accurate relation. 

In this research, relation between wave amplitude A and 
wave constant G was formulated without working on the 
assumption of small and long wave and without 
eliminating wave constant G. As a result, an equation 
where wave constant G together with wave amplitude A 
as its variables was formulated. The relation was 
formulated using KFSBC together with momentum 
equation. 

To simplify the calculation, the equation was developed 
by maintaining an assumption that the wave is sinusoidal, 
where wave height H is twice wave amplitude or H = 2A. 
The existing understanding is that sinusoidal wave is only 


for small amplitude wave. This research experienced no 
obstacle as a result of the working on sinusoidal wave 
assumption and quite accurate calculation result was 
obtained. 


II. VELOCITY POTENTIAL FOR SLOPING 
BOTTOM 

hi this research, velocity potential at sloping bottomfrom 
Hutahaean (2008) was used with initial form as follows 
0 = Ge kh p(z)coskxsinat 

G is wave constant, k is wave number, a is angular 
frequency, a = ^ , T is wave period, z is elevation at 
vertical axis with z = 0 at still water level surface. In this 
research, a new constant G = Ge kh was formed and 
bearing in mind the separation between G and e kh did not 
provide any benefit at all; therefore, to make it more 
concise,incorporation was done, then the potential 
velocity equation becomes: 

0 = Gfl(,z) coskxsinat .(1) 

/?(z) = ae kCh+z) + e ~ kih+z) .(2) 

^(z) = ae k(h+z) - e ~ k(h+z) .(3) 

/ 1+ 3£ i_3£\ 

~ M a* . <ix \ 

a = - 2 \-^ + ffB) .< 4 > 

\ a* a*/ 

— is bottom slope. At sloping bottom there will be — and 
jj—. Related to the changes in the value of the wave 
constant, there are two conservation laws in the velocity 
potential equation, i.e. wave number conservation and 
energy conservation, 
a. Wave Number Conservation 

At the execution of Laplace equation with variable 
separation method, as it is with coshk (h + z ) at velocity 
potential for horizontal bottom, it has been determined 
that /? (z) is just a function of z, hence: 

2^=0 .(5) 

ax 

= 0 .(6) 

ax 

Atz = 0, applies 

f- = 0 .(7) 

or 

a k k ah 

= . .(o) 

ax n ax 

For z = t) where t] is the elevation of water wave surface 
against still water level. 
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3% 


■= 0 


At the characteristic point where coskx = sinkx = 
cosat = sinat , hence n = - , where A is wave 

amplitude. This will be explained in other part. Wave 
number conservation equation becomes 
ak( h+ f) 


0 


•(9) 


i 3 k k 3A , a h _ .. „. 

■ +- +k— = 0 ( 10 ) 

2 ^ a* 2 a% a% 


(h+±y 

V 2/; 

Equation (10) can be called continuity equation or mass 
conservation equation at waves experiencing changes in 
water depth, wave number k and wave amplitude A. This 
equation cannot stand by itself, it has to meet certain 
limitation conditions, in this case Kinematic Free Surface 
Boundary Condition (KFSBC). 

b. Energy conservation 

Other characteristic contained in (1) is energy 
conservation that will be fonnulated as follows. 

The velocity of water particle in horizontal u direction at 
the direction of axis-x is, 

u = — — = GkB (z) sinkxsinat — 

3x 

— B (z) coskxsinat .(11) 

a* r 

Differential of this particle velocity equation in horizontal 
direction against axis-x, where there are changes in G and 
k is 
3u 

— = Gk 2 B(.z) coskxsinat 
3x 

3 Gk 

H- B(z) sinkxsinat 

3x 

+ — kB(z) sinkxsinat — -—— Biz) coskxsinat 
Particle velocity in vertical direction is, 

w =-= —GkBi(z) coskxsinat 


a z 


..( 12 ) 


The differential of this equation against vertical-z axis is 
3 w 

— = —Gk 2 Biz) coskxsinat 
3 z 

_ _ , . 3u 

Mass conservation law or continuity equation is-1- 

a% 

— = 0. Substitute equations —and — , and work on the 

a z 3x a z 

characteristic point , i.e. a point where coskx = sinkx = 

cosat = sinat, 

at ac a 2 c 

G -h 2 — k - - = 0 

a% a% a% 2 


,.(13) 


The equation can be written as, 


a ‘■g 


■■(14) 


r ak a g 

, G -h 2 —k 

3x 2 3x 3x 

Either (13) or (14) can be called mass conservation 
equation, but considering that at G there is energy 


dimension, then it can be called energy conservation 
equation. At either long wave or small amplitude an 

g 2 G 

assumption of —— = 0 can be done, so the following 
a% z 

relation was obtained, 

.(15) 


ak 

^ = -&-G 

3x 2k 

or 


f = -hGfc .(16) 


^ .( 17 > 

3k 


..(18) 


at (15), (16) and (17) can be approached with ( 8 ) i.e. 

3x 

— — Therefore, in this case energy conservation 

ax h ax 

equation is connecting changes in G with changes in 
water depth h and wave number k. In this research energy 
conservation equation with higher degree of accuracy was 

g3 q g3 q 

used, i.e. = 0. Equation for —was obtained by 
differentiating (14) against horizontal-x axis and at the 

g2 q g2 ^ 

element —7 substitute with (14), and is ignored since 

it is very small, a fonn of energy conservation like (16) 

was obtained, i.e. — = — ufcG with different coefficient 
ax r 

of change p. 

2— 

With relation (16), hence particle velocity at the direction 
of horizontal-x axis becomes 
u = Gk (sinkx + \icoskx) B(z)sinat....( 19) 

III. THE METHOD OF WAVELENGTH 
CALCULATION 

The first equation connecting G and k is KFSBC: 

Y — = w„ ~ u„ — ....(20) 

'at 71 71 ax 

Y is weighting coefficient at weighted total acceleration 
equation with a value of 2.784-3.160 (Hutahaean (2019a)) 
and 2.483 for Hi/ and 2.202 for Hi/ in Hutahaean 

(2019b). This research used the value ofy = 2.483. 
Substitute (12) and (19) that was done at z = r ]to (20), 
where the right side is multiplied with -, 
a?7 Gk 

at ya 

3m 


/ 

-a (&(? 7 ) 


coskxsinat 


+ 


dr]\ 

(sinkx + \icoskx) B (r]) sinat —) 

3x/ 


-{^BiW) + (l + j acoskxsinat 


At the characteristic point, this equation can be written as 
an Gk ( 

Tt = ~^ 

Waveanphtude A is defined as, 

A = ^(/? 1 (t 7 ) + (l + p)/?(j 7 )g).( 21 ) 

— = —Ao coskxsinat ....(22) 
at v 
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rj = Acoskxcoscrt ....(23) 

— = — kAsinkxcos at ....(24) 

dx 

By working on the assumption at (22), (23) and (24), 
there is an assumption that water wave surface equation is 
sinusoidal, where wave height H = 2 A. From (23), at the 

A 

characteristic point rj = -, then (21) becomes 

' 4 = ^( ft (;)- (1 + ^© T >. (25) 

This equation is the first equation for G and k calculation 
with wave amplitude A as input, 


...(26) 

As the second equation is surface momentum equation 

where convective velocity is ignored. 

3u v a., 

= . (27) 

From (19) 

u ); = Gkisinkx + \icoskx) flip) s inert 

3u, ; 

y-^- = Gkaisinkx + \icoskx')/3(j]')cos<7t 

At the characteristic point the second equation was 
obtained, i.e. 

f 2 (G, fc) = yoG (l + p)/? (-) - gA 


.(27) 

(27) can be formed into equation for Gand it is substituted 
to (26), so an equation is formed just for wave number 
k. In this research, the two equations are done 
simultaneously, so the values of G and k are obtained 
which meet KFSBC and momentum equation, with input 
a, h and wave amplitude A. The calculation was done 
using Newton-Rhapson method, where this iteration 
method needs initial price of iteration. As the initial price, 
for wave number k equation from Hutahaean (2019a) was 
used, i.e. 

Y 2 ° 2 = gk{ 1 -(f)) .(28) 


This equation is a quadratic equation of wave number 
fcthat can be completed using simple method to find the 
root of a quadratic equation. After the value of wave 
number was obtained, the initial estimation value of Gcan 
be calculated using (25) where wave amplitude as input, 

c Jim . <29) 

3.1. The result of wavelength /.calculation 
Wave number k can be changed into wavelength with a 
simple relation, i.e. L = f. The result of wave length 
calculation for various wave periods in a number of water 
depths is shown on Fig.l, where the calculation was done 
at bottom slope— = —0.005, weighting coefficient y = 

3% 

2.483 , and wave amplitude A = 0.6. m 



— T = 8 sec - T = 9 sec -T = 10 sec 

Fig.l: Graph of wave length against water depth 

Fig.l shows that wave with wave period of 8 sec, 
wavelength constant starts at water depth 9.0 mor deep 
water h 0 = 8 m with L 0 = 14.0 m, whereas at wave 
period 9 sec., h 0 = 10.0 m with L 0 = 18.4 m, at wave 
period 10 sec. h 0 = 12 m, with L 0 = 23.2 m The wave 
length is quite realistic, quite in line with what exists in 
the nature. 

At the equations for calculating wave number, there are 
wave amplitude as the variables, hence the wave length 
that was produced will be determined by the size of the 
wave amplitude A. Fig.2 shows graph of wavelength 
against water depth using wave with wave period of 8 
sec., wave amplitude A, 0.60 m, 0.80 m and 1.0 m It is 
shown that the bigger the wave amplitude, the shorter the 
wave length. The calculation was done using — = 
— 0.005, weighting coefficient y = 2.483. 


15 
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- A = 0.60 m - A = 0.80 m -A = 1.0 m 

Fig.2: The influence of wave amplitude on wave length 

With the presence of the influence of wave amplitude on 
wave length, the calculation of wave length change from a 
water depth to a shallower water depth should have been 
calculated along with shoaling analysis. There is also the 
influence of bottom slope on wave length which will be 
discussed on the next research due to space limitation. 
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Fig. 3: The value of wave constant G 


Fig.3 shows the value of wave constant G for wave with 
wave period of 8 sec., where the influences of wave 
amplitude and water depth are visible. At the deep water, 
the influence of wave amplitude is not that big where the 
value of G at the three wave amplitude values is almost 
the same. The difference is visible at shallow water, 
where the bigger the wave amplitude, the bigger the value 
of G but at big wave amplitude, i.e. 0.80 m and 1.0 m 
where the difference is very small. Considering that 
basically particle velocity is a function of wave constant 
G, then if wave constant G is eliminated with wave 
amplitude, then a very accurate relation between wave 
constant G and wave amplitude is needed in order to 
obtain an accurate particle velocity. 


IV. SHOALING AND BREAKING MODEL 


In the journey to shallower water, there are 3 (three) wave 
parameters with changing values, i.e. wave number k, 
wave constant G and wave amplitude A, where in this 
case wave amplitude A has been absorbed as variable of a 
wave. Therefore, at shoaling and breaking analysis there 

are 3 (three) unknowns, i.e. —, — dan — where the 
' ax ax ax 

changes occurred as a result of changes in water depth h. 
There is an unchanged wave characteristic, i.e. 



Index 0 shows the value at deep water. The conservation 
of the two characteristics is the consequence of wave 


number conservation law, where 


^W = 0and 5MW = 

3x 3x 


o. 


. r . . 3 k 3 G . 

As a governing equation of changing equation —, — and 

3x cix 

Tj^-are wave number conservation equation, i.e. (10) 

[h+-) — + - +k— = 0 .(10) 

v 2 / ax 2 ax ax 

The next equation is KFSBC equation in the form of an 
equation for wave amplitude A (25), 

“). <25 » 

Defined, 


/= 0*1(7)-Cl+ i0/?(f)“) 

•(31) 


..(30) 


Gk 


yo 


f 


Wave amplitude equation (31) was differentiated against 

horizontal-x axis 

aA (ac , „ ak\ r 

i"’5T = t' I + c ;r )f .< 32 > 

3^4 

Then—is substituted to (10) 

3x 


/= 0 


( ^ 

,a k 

a h 

k 

/3G 

a k\ 

(h +- 

I h k h 

- 1 

— k + G — 

V 2) 

dX 

dX 

2y<j 

\dX 

dX/ 


..(33) 


Hence (33) meet KFSBC. For (27) to meet energy 
conservation equation, then — is substituted to (16), — = 

3x 3x 

— \ikG, to obtain changing equation of wave number — 

3% 

that meets wave number conservation equations,, KFSBC 
and energy conservation. 

( h+ ±) + cM-Y± = _ k ajL + H£*!£ . (34) 

' 2 ' 2vcr / 3x 3x 2 yo 


At (34), there is variable \i which is a function of 


3 k 


. 3 k 


Therefore, the calculation of —with (34) was done using 

3x 

iteration, i.e. the first step, —was calculated with (8),— = 

3x 3x 

— the value of uwas calculated with (18). Then, 
h ax ^ v ' 

—was calculated with (34) andpwas recalculated with 

3x 

(18). This step is repeated over and over again until a 
stable value of — where generally 5-6 iterations have 
obtained stable value of — . After a stable value of— was 

3x 3x 

obtained, — was calculated with (16) and — with (10). 

3x 3x 


h x + Sx —where—is 
x ax ax 


Then, the value of variable at the point x = x + bxwas 

calculated at the depth of h x+Sx 

negative, using Taylor series. 

3 k 

kx+8x ~ k x T Sx — 


J x+Sx 


^x+Sx 


Gy SX ' 


: A x + Sx ■ 


3x 

3G 

3x 

3A 

3x 


The calculation was done from the deep water until 
coastal water, until breaker point and afterward 
4.1. Breaking Characteristics 
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At the equation for the wave amplitude, i.e. (30) and (31), 
there is breaking condition when / = 0, 

(/S 1 (^)-(1 + m)/?(^)“) = 0....(35) 

fi 1 (y) and /? (y) are values at deep water, where 

/?j = /? ("^), so that if (35) is divided by /? (-^)the 


result is 


.(36) 


^b ^b _ r 

2 (l + n) 

In (36), there is the influence of bottom slope on the 
breaking parameter that was absorbed at the value of p. 
The formulation process has stated that wave is 
sinusoidal, i.e. in (22). (23) and (24), then relation A b = 
applies, therefore breaker length index is obtained. 

2 
Lb 


Hb 

2 

Hb 


(l+ [J.) 7T 


.(37) 


In (37), there is the influence of bottom slope on breaker 
height, whereas the previous section has shown that there 
is the influence of bottom slope on wavelength. So, in 
general, bottom slope will have an influence on the 
breaking wave. However, more detail discussion will be 
done in the next research. 


4.2. The result of Shoaling-Breaking Model 

As an example, a model was executed for wave with 

wave period T = 8 sec., deep water wave height is H 0 = 

2.32 m, bottom slope — = —0.005, weighting 

3x 

coefficient y = 2.483. As seen in Fig.4, breaking occurs 
at breaker depth h b = 3.66 m, with breaker height H b = 
2.95 m. 



Fig. 4: Shoaling-breaking, wave period T = 8 sec.. H 0 = 
2.30 m. 


The comparison with breaker height index equations is 
presented on Table (1). As comparator, the average values 
of 5 empirical breaker index equations that were obtained 
from laboratories experiments were used, i.e. Komarand 
Gaughan (1972), Larson, M. and Kraus, N.C. (1989), 
Smith and Kraus (1990), Gourlay (1992) andRattana 


Pitikonand Shibayama (2000). Whereas comparator for 
breaker depth, equation from SPM (1984) was used. 

Table. 1: Comparison of the result of the model with 


empirical equation 


T 

W 0 

H b 

(m) 

h b (m) 

(sec) 

(m) 

Model 

BHI 

Model 

SPM 

8 

2,32 

2,95 

2,71 

3,66 

3,38 

9 

2,94 

3,73 

3,43 

4,63 

4,28 

10 

3,63 

4,6 

4,23 

5,71 

5,28 

11 

4,39 

5,57 

5,12 

6,91 

6,39 

12 

5,22 

6,63 

6,09 

8,22 

7,61 

13 

6,13 

7,78 

7,15 

9,65 

8,93 

14 

7,11 

9,02 

8,29 

11,19 

10,36 

15 

8,16 

10,36 

9,52 

12,85 

11,89 

16 

9,29 

11,79 

10,83 

14,62 

13,53 


Note: BHI: average from 5 (five) Breaker Height Index 
equation 


The comparison between the result of the model with 
empirical breaker height equation is shown on Table (1). 
As deep water wave height H 0 , 0.9 x H 0 _ max was used in 
every wave period, where H 0 _ max= —p—; m (Hutahaean 

(2019b)), whereas deep water depth h 0 = , bottom 

k 0 

slope — = —0.005, weighting coefficient y = 2.483.. 
Breaker height H b from the model is bigger than breaker 
height from BHI, with an increasing pattern of differences 
with the increase in wave period.Similar differences 
pattern occur at breaker depth h b , but with not so big 
differences. For more clear information see Fig .5 for 
breaker height comparison and Fig.6. for breaker depth 
comparison. 



-- Model - BHI 

Fig.5: Comparison of breaker height 
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— model -SPM-1984 

Fig. 6: Comparison of breaker depth 


Table.2: Comparison of breaker height index — 

h b 


T 

(sec) 

h„ 

Model 

Model 

BHI 

Lb 

(m) 

Hb 

Lb 

8 

0,81 

0,8 

4,6 

0,64 

9 

0,81 

0,8 

5,82 

0,64 

10 

0,81 

0,8 

7,19 

0,64 

11 

0,81 

0,8 

8,7 

0,64 

12 

0,81 

0,8 

10,35 

0,64 

13 

0,81 

0,8 

12,15 

0,64 

14 

0,81 

0,8 

14,09 

0,64 

15 

0,81 

0,8 

16,17 

0,64 

16 

0,81 

0,8 

18,4 

0,64 


—^comparison is shown in Table (2), which shows that the 
hb 

result of the model is quite close with the result of 

empirical equation. Breaker steepness cannot be 

compared since there were wave-length differences, 

where breaker length model L b is quite short, quite in 

accordance with what exist in the nature. Breaker 

steepness —is quite in accordance with analytical 

Lb 

equation (37), i.e. — « -. 

Lb n 

V. CONCLUSION 

Wave dynamic calculation using wave constant G can be 
done easily and provide a quite good calculation result. 
The wave constant execution enables energy conservation 
equation execution which is a relation between wave 
constant changes with wave number changes. 

The execution of sinusoidal wave assumption, simplify 
the correlation of calculation result in the form of wave 
amplitude with wave height, where generally the 


ISSN: 2349-649S(P) / 2456-1908(0) 

information needed is wave height, hi addition, the 
execution of sinusoidal wave assumption in this research 
found no difficulty in the use of even big wave amplitude. 
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